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ABSTRACT 

The Kepler mission is dramatically increasing the number of planets known in multi-planetary systems. Some 
pairs of close adjacent planets have orbital period ratios close to resonant values, with a tendency to be larger 
than required for exact first-order mean-motion resonances . This intriguing feature has been shown to be a 
natural outcome of energy dissipation in resonant planetary pairs with conserved angular momentum, and tidal 
circularization due to star-planet interactions has been highlighted as a possible mechanism for such energy 
dissipation. However, many multi-planetary systems have orbital periods larger than ten days, for which tidal 
circularization is unlikely to provide efficient energy dissipation. Gravitational interactions between partial 
gap-opening planets and their parent protoplanetary disk may instead provide efficient dissipation. We show 
that disk-planet interactions can lead to resonant planet pairs with period ratios in significant excess to that of 
first-order resonant period ratios. It could help account for the diversity of period ratios in Kepler's planet pairs 
comprising super-Earth planets with orbital periods greater than about ten days. Fast disk-driven migration of 
sub-jovian planets may alternatively form planet pairs evolving in a common gap, where the planets' convergent 
migration may stall with a variety of period ratios. This mechanism is found to reproduce the period ratio of 
Kepler-46b/c. 

Keywords: accretion, accretion disks — hydrodynamics — methods: numerical — planetary systems: forma- 
tion — planetary systems: protoplanetary disks 



1. INTRODUCTION 

Planetary astrophysics is undergoing an epoch of explosive 
growth driven by the discovery of about 850 exoplanets in two 
decades. Continuous improvement in detection techniques 
will uncover many more planets in the near future, which 
are too small or too far from their host star to be detected at 
present. Observations suggest that planet formation is ubiqui- 
tous, with as many as 80% of nearby Sun-like stars harboring 
one or more Earth- to Neptune-analogues on close-in orbits 
dMayor et al.ll2009h . These prolific discoveries have revealed 
the amazing diversity of exoplanetary systems, and indicated 
that planets may form under a wide range of conditions and 
have considerable mobility. 

The diversity of exoplanets has been particularly high- 
lighted by the Kepler's space missi on. To date, Kepler h as de- 
tected over 2300 planet candidates (Batalha et al. 2012), from 
which 105 are confirmed planets (http://kepler.nasa.govi. 
About one third of Kepler's candidates are associated with 
multiple transiting systems (765), and the vast majority 
of them are expecte d to be real multi-planetary systems 
dLissauer et al.l 120121) . Some of these systems show re- 
markable properties, like the five coplanar, alternating rocky 
and icy planets or biting Kepler-20 in less than 80 days 
dGautier et al.l 120121) . Interactions between planets and their 
parent protoplanetary disk are likely to have played a major 
role in shaping such compact system. Another striking fea- 
ture of Kepler's multi-planetary systems is the large number 
of planet pai rs that are not in me an-motion resonance (here- 
after MMR, iLissauer et al.ll201ll see also Figure [T). How- 
ever, planet pairs show some tendency to have period ratios 
slightly greater than resonant values. For instance, there are 
almost twice as many planet pairs with period ratios between 
2.0 and 2.1 than between 1.9 to 2.0. This intriguing feature 
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has been shown to naturally arise from tid al orbital circular- 
ization of close-in resonant planet pair s ( Papal oizoul 120 1 It 
iLithwick & Wul l20ll iBatvgin & Morbidellil 12013k How- 
ever, as emphasized in FigureQ] in many multiple-planet sys- 
tems the inner planet has an orbital period exceeding ten days. 
For these systems, tidal circularization is unlikely to have 
caused significant divergent evolution of the planets orbits. 

An example of a planetary system with a period ratio 
slightly greater than resonant is the Kepler-46 system (pre- 
viously known as KOI-872). The 0.9 Solar-mass star is sur- 
rounded by a 0.8/?j transiting planet on a 33.6 day period 
(Kepler-46b). Large transit time variations (TTV) of Kepler- 
46b are caused by the presence of an outer planet, Kepler-46c, 
whose mass (0 .37M|) and orbital peri od (57.0 days) have been 
determined bv lNesvornv et ail (120121) . The same authors have 
also reported the detection of a 1 .77?® planet candidate on a 
6.8 day period (Kepler-46d, not yet confirmed). The main 
properties of the Kepler-46 planetary system are summarized 
in Table [T] The low-eccentricity and quasi-coplanar orbits in 
this system suggest that the planets reached their present loca- 
tion most probably through disk-planet interactions. Still, the 
ratio or orbital periods between Kepler-46c and Kepler-46b 
is w 1.696, which is slightly greater than 5/3. Capture into 
the 5:3 MMR is not a natural outcome of planet migration 
driven by disk-planet interactions. It requires a rate of conver- 
gent migration that is large enough to cross the 2: 1 MMR, but 
not too large so that the planets do not reach their 3:2 MMR. 
Capture into a second-order resonance also requires finite ec- 
centricities, which disk-planet interactions tend to damp quite 
efficiently. Should disk-planet interactions have led Kepler- 
46b and Kepler-46c to reach their 5:3 MMR, it is unclear how 
their period ratio would have then increased to their present 
value. 

This paper is organized as follows. In Section|2] we present 
results of hydrodynamical simulations modeling the early 
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Figure 1. Ratio of orbital periods of planet pairs among all candidate multi- 
planetary systems detected by Kepler, shown as a function of the orbital pe- 
riod of the inner planet (in days, y-axis). The circles' colour is varied accord- 
ing to the physical radius of the inner planet in each pair, expressed in Jupiter 
radii. Vertical dashed lines show the location of the 2:1, 5:3, 3:2 and 4:3 
mean-motion resonances. The top panel is an histogram of the orbital period 
ratio. Data were extracted from http://planetquest.jpl.nasa.gov/kepler/ 



Table 1 

Kepler-46 planetary system jNesvorny et a l.120121) 



Mass or radius Period [days] Eccentricity 



Kepler-46 


0.9M : , 






Kepler-46<f 


1.7R e 


6.8 


(assumed) 


Kepler-46b 


« 0.8R, 


33.6 


<0.02 


Kepler-46c^ 


0.37Mj 


57.0 


0.015 



"Unconfirmed 

"Inferred from TTV signal of Kepler-46b. 



evolution of Kepler-46b and Kepler-46c as they were embed- 
ded in their parent protoplanetary disk. These Saturn-mass 
planets are expected to open a partial gap around their orbit. 
A rapid convergent migration causes the planets to merge their 
gap and to evolve into a common gap. The planets' conver- 
gent migration is found to stall with period ratios between 
1.5 and 1.7, depending on the disk's density profile inside the 
common gap. The observed period ratio can be reproduced 
without the planets being in resonance. Furthermore, we show 
that disk-driven migration of partial gap-opening planets may 
lead to significant divergent evolution of resonant planet pairs 
over typical disk lifetimes. This mechanism is illustrated for 
super-Earth-mass planets in Section|3]with both hydrodynam- 
ical simulations and customized three-body integrations. Di- 
vergent evolution of resonant gap-opening planet pairs medi- 
ated by disk-planet interactions could partly account for the 
diversity of period ratios in Kepler's multi-planetary systems. 



2. HYDRODYNAMICAL SIMULATIONS OF THE 
EARLY EVOLUTION OF THE KEPLER-46 
PLANETARY SYSTEM 

We model in this section the early evolution of the Kepler- 
46 planetary system when the planets are embedded in their 
parent protoplanetary disk. Our study does not address the 
formation of the planets. We focus instead on the planets or- 
bital evolution driven by disk-planet and planet-planet inter- 
actions. For this purpose, two-dimensional hydrodynamical 
simulations of plan et-disk interac tions were carried out using 
the code FARGO (lMasselll2000l) . The physical model and 
numerical setup o f the hydrodynamical simulations are de- 
scribed in Section 12.11 Results of simulations follow in Sec- 
tion O 

2.1. Physical model and numerical setup 

Disk model — We adopt a two-dimensional disk model for 
the protoplanetary gaseous disk in which the planets in the 
Kepler-46 system formed. Disk-planet interactions determine 
the rate of convergent migration of the planets as well as the 
damping rate of their eccentricity. The migration and eccen- 
tricity damping rates are sensitive to the disk's surface den- 
sity, temperature, turbulent viscosity and radiative properties. 
These make the parameter space particularly large, and the 
need to assess the planets' resonant dynamics over typically 
a thousand planet orbits has led us to make simplifying as- 
sumptions for the disk model and to fix some of the disk pa- 
rameters. First, gas self-gravity is discarded for simplicity as 
the disk masses adopted throughout this study are low (the 
Toomre-Q parameter associated with the unperturbed density 
profile does not go below 7). Also, a locally isothermal equa- 
tion of state is used where the vertically-integrated gas pres- 
sure P and density E satisfy P = Ec^, with the gas sound 
speed c s being specified as a fixed function of the cylindri- 
cal radius r. The sound speed is related to the disk's pressure 
scale height H through H = c s /£2k with £2k the local Keple- 
rian angular velocity. The (fixed) temperature profile is taken 
proportional to r , which leads to a uniform disk's aspect 
ratio h = H/r. We take the generic value h = 0.05. The ef- 
fect of turbulence is modeled by a constant shear kinematic 
viscosity v. We take v = 1.1 x 10~ 5 in standard code units 
(defined below) which translates into a viscous alpha param- 
eter a — 4.5 x 10~ 3 at the initial location of the inner planet. 
This value of a is typical of the magnetically active regions of 
protoplanetary disks with magneto-hydrodynamic turbulence 
driven by the magneto-rotational instability. The present or- 
bital periods of the planets in the Kepler-46 systems imply that 
the planets have probably interacted with such active regions 
in their parent disk. Note however that the disk's magnetic 
field is not included in our simulations. We also point out 
that the planets we will consider are of large enough mass so 
that small-scale turbulent fluctuations can be safely ignored 
when modeling the interaction with their parent disk. In other 
words, the viscous diffusion approximation is expected to be a 
reasonabl e approximation for the large planet masses we w ill 
consider dNelson & Papaloizoull2003t iBaruteau et al.ll20Tl1) . 

The disk is set up in radial equilibrium, with the centrifugal 
acceleration and the radial acceleration related to the pressure 
gradient balancing the gravitational acceleration due to 
the central star. The initial gas surface density is taken 
proportional to r~ 1//2 . Its value at the initial location of the in- 
ner planet is taken as a free parameter, which we denote by Eq. 
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Planets parameters — The planets are assumed to have 
already formed at the beginning of our simulations. They 
therefore start with their observed mass. We only model the 
orbital evolution of the outer two planets in the Kepler-46 
system: Kepler-46b (inner planet) and Kepler-46c (outer 
planet). The 1.7/?© super-Earth at 6.8 days (Kepler-46d), if 
confirmed, is assumed to have had a negligible impact on the 
evolution of Kepler-46b and Kepler-46c. Kep ler-46b has a 
physical radius « 0.8/?j (Ne svornv et al.l 120121) . The mass- 
radius diagram of exoplanets detected by transiting methods 
indicates that the mass of Kepler-46b most probably lies be- 
tween 0.2Mj and 0.6Mj. Its planet-to-primary mass ratio (q{) 
is therefore varied in the range [2.2 x l(T 4 -6.6x 1(T 4 ]. The 
TTV signal of Kepler-46b constrains t he mass of Kepler-46c 
to be « 0.37Mj (N esvornv et al.ll2012l) . Its planet-to-primary 
mass ratio is thus fixed to q = 4 x 10~ 4 . The initial orbital 
radius of Kepler-46b, rj, is taken to be the code's unit of 
length. That of Kepler-46c is arbitrarily set to r Q = 1.8n, 
so that the ratio of orbital periods m 2.4 initially. For the 
disk aspect ratio and viscosity taken in our study, the two 
planets are expected to open a partial gap around their orbit 
dLin & Papaloizou|[T993l:ICrida et al.ll2006l) . 

Numerical setup — The hydrodynamical equations are solved 
in a cylindrical coordinate system {r, <p} centred on to the star 
with r G [0.2 — 2.6] r\ and <p e [0,2tt]. The computational grid 
has /V r = 400 zones evenly spaced along the radial direction 
and N s = 800 azimuthal sectors. The frame rotates with 
the Keplerian frequency at the inner planet's location, and 
the indirect terms that account for the acceleration of the 
central star by the disk and the planets are included in the 
equations of motion. A standard outflow boundary condition 
is taken at the grid's inner edge, while damping is used in a 
so-called wave killing-zone extending from r = 2.2 to r = 2.6 
(disk quantities being damped towards their instantaneous 
axisymmetric profile). To avoid a violent relaxation of the 
disk due to the sudden introduction of the planets, the mass 
of the planets is gradually increased over 10 orbital periods. 
The gravitational potential of the planets is smoothed over a 
softening length, e, equal to 0.6// with H evaluated at the 
planets position. The calculation of the force exerted by 
the disk on the planets excludes the content of the planets' 
circumplanetary disk , the size of which is about 60% of the 
planets' Hill radius dCrida et al.l 120091) . The Hill radius of 
each planet is resolved by about 10 grid cells along each 
direction initially. 

Code units — Results of simulations are expressed in the stan- 
dard units of disk-planet calculations: the mass unit is the 
mass of the central star (denoted by M*, and equal to 0.9M© 
for the Kepler-46 system). The length unit is the initial orbital 
radius of the inner planet (r0, and the time unit is the initial 
orbital period of the inner planet (r or b) divided by 2n. When- 
ever time is expressed in orbits, it refers to the orbital period 
at the initial location of the inner planet. 

2.2. Results of hydrodynamical simulations 

The free parameters in our hydrodynamical simulations are 
the mass of the inner planet and the unperturbed disk's den- 
sity profile, which control the rate of convergent migration of 
the planets and their eccentricity. We consider three different 
masses for the inner planet: M mneY = 0.2Mj, 0.4Mj and 0.6Mj 
(recall that the mass of the outer planet is fixed to 0.37Mj). 



For each planet mass, a series of 9 simulations was carried 
out varying T.q, the unperturbed surface density at the ini- 
tial separation of the inner planet. We took £o in the range 
[2—10] x 10~ 4 , which translates into an initial disk mass vary- 
ing from 3 .5 x 1( T 3 M* t o 1.8 x 10~ 2 M*. Our results are de- 
scribed in Section 12.2.1 1 To obtain evolution of the planet's 
orbital elements to a steady state, a simple model for di sk dis- 
persal is included and its results presented in Section 12.2.21 
The robustness of our results to varying som e phys ical and 
numeri cal par ameters is discussed in Section 12.2.31 Lastly, 
Section l2.2.4l shows that two planets that do not directly inter- 
act gravitationally can still exhibit resonant coupling through 
the indirect mediation of the disk. 

2.2. 1 . A few illustrative cases 

Figure [2] shows the typical outcomes of our hydrodynami- 
cal simulations for the Kepler-46 system. From left to right in 
the figure, the panels display the time evolution of the plan- 
ets' semi-major axes, orbital period ratio, eccentricities and 
resonant angles (Aj and <»j. denote the planets' mean longi- 
tude and longitude of periastron indicated, respectively, and 
the i and o subscripts refer to the inner and outer planet). The 
first row from the top shows results with M mner = 0.4Mj and 
Eo = 3 x 10 4 . Due to the opening of a partial gap around 
each planet's orbit, the relative convergent migration of the 
two planets does not occur at a constant rate. In particular, the 
outermost planet goes through a stage of fast type Ill-like mi- 
gration before slowing down upon approaching the 2: 1 MMR 
with the inner planet. The planets maintain their 2:1 MMR 
from 400 orbits onwards, as can be seen from the libration 
of the two resonant angles associated with the planets' 2:1 
MMR. The resonance capture increases the inner planet's ec- 
centricity to about 0.03. The increase in the outer planet's 
eccentricity is much more modest, reaching a few xl0~ 3 . 
Damping of both eccentricities follows from the gravitational 
interaction with the background gas disk. Interestingly, we 
see that the ratio of orbital periods rises significantly after 
the resonant capture: it increases from 2.02 up to 2.15 in 
about 1500 orbits. Note that resonant coupling is maintained 
throughout this evolution. 

The increase in the orbital period ratio is a consequence 
of the planets relative migration being nearly stalled while 
their eccentricity is continuously damped by the disk. It can 
be sketched by saying the two-planet system evolves at ap- 
proximately constant total angular momentum with total en- 
ergy dissipation provided by disk-driven eccentricity damp- 
ing. Alternatively one can consider the system to be evolv- 
ing with a slowly decreasing ratio of circularization time to 
migration time. It is known that if this ratio is fixed, a reso- 
nant planet pair attains equilibrium eccentricities that are pro- 
portional to the square root of this rati o, after which the sys- 
tem evolves in a self-similar manner (|Nelson & Papaloizoul 
l2l)02t IPapal oizou & Szuszkiewiczll2005l) . Suppose now this 
ratio slowly decreases to zero on account of the migra- 
tion time becoming arbitrarily large while the circularization 
time remains fixed. Then the eccentricities will slowly de- 
crease as the system moves through a sequence of quasi- 
steady states. Moreover, as the eccentricities are resonantly 
forced, they are inversely proportional to the deviation of 
the ratio of the semi-major axes from the strictly resonant 
value. Accordingly, assuming that the period ratio starts 
slightly larger than the strictly resonant value, the period ra- 
tio must steadily increase with time. If the migration rate 
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Figure 2. From left to right: time evolution of the planets' semi-major axis, orbital period ratio, eccentricity and of the resonant angles relevant to the planets' 
period ratio. Each row corresponds to a different hydrodynamical simulation of the Kepler-46 system. The mass of the inner planet (A/j nnel ) and the disk s 
unperturbed surface density at r = 1 (Eo) are shown in the panels' top-right comer. The location of some mean-motion resonances are displayed by dashed curves 
in the semi-major axes panels. The first row shows an example of slow convergent migration towards the 2:1 mean-motion resonance followed by rapid resonant 
repulsion. The second and third rows show cases where convergent migration stalls with orbital period ratios near 1.7 (observed value) and 1.6. The lower row 
shows a case where the planets get close to the 3:2 mean-motion resonance followed by slow resonant repulsion. 



is reduced to arbitrarily small values, the period ratio can 
increase to large values, and the system behaves like one 
that is undergoing circularization, and therefore losing or- 
bital energy at fixed angular momentum. The energy is dis- 
sipated through interaction with the disk. This mechanism 
has been recently studied in the context of the late evolu- 
tion of planetary systems where energy dissipation arises 
from orbital circularization due to tidal interaction with the 
central star (iPapaloizou & Terquemll2010t lPapaloizoull20TTt 



iLithwick & WulEOia iBatvgin & Morbidellill2013l) . We re- 
fer to the mechanism discussed in this paper as disk-driven 
resonant repulsion, fol lowing the terminology introduced by 
ILithwick & Wul (120121) in the context of tides. 

No or very little resonant repulsion occurs in the other pan- 
els in Figure [2] The second row from the top displays results 
for Mi nner = 0.6M] and Eo = 6 x 10 4 . This case is particu- 
larly interesting as it shows that convergent migration may 
stall with an orbital period ratio very close to the observed 
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Figure 3. Disk's surface density obtained at 400 orbits for M- mn 
and Eo = 3 x 10~ 4 (upper panel), and for Mj mlcr = 0.6Mj and £o : 
(lower panel). 
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TTV value (« 1.7). As the planets cross the 2:1 MMR at 
about 100 orbits, their eccentricities increase and convergent 
migration slows down. The decreased rate of convergent mi- 
gration is not a consequence of resonant repulsion (the 2:1 
resonant angles do not librate) but is rather due to the planets 
carving a common gap in the background disk. This is illus- 
trated by the two panels in Figure[3] which compare the disk's 
surface density at 400 orbits for the two aforementioned sim- 
ulations: Mi nn er = 0.4Mj and Eo = 3 x 10 4 (top panel), and 
Mnner = 0.6M] and Eo = 6 x 10 4 . While in the former case 
each planet opens a gap with a narrow ring of gas left in be- 
tween, the latter case shows that both planets end up in a com- 
mon gap. Back to the second row of panels in Figure [2] we 
see that the orbital period ratio remains approximately sta- 
tionary from about 400 orbits, and so do the eccentricities. 
The time-averaged eccentricity of the inner planet is close to 
0.01, and that of the outer planet is 0.025. These values are 
consistent, albeit slight larger than those inferred from obser- 
vations (see Table [!}. The third row of panels in Figure |2]is 
obtained for the unperturbed disk's density slightly increased 
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Figure 4. Summary of the results of hydrodynamical simulations for Kepler- 
46: ratio of orbital periods versus the unperturbed disk's surface density at 
the initial location of the inner planet (Eo). Results are shown for the three 
masses assumed for the inner planet. 



to Eo = 7 x 10 4 . The planets' orbital evolution is very similar 
to the previous case with Eo = 6 x 10 4 , except that planets get 
captured in their 3:2 MMR with an orbital period ratio « 1 .6, 
which can happen because of the very small eccentricities in 
this case. 

The lowest row of panels in Figure [2] shows another out- 
come of the simulations (M; nner = 0.2Mj and Eo = 8 x 10 4 ). 
Here, convergent migration is rapid enough that the planets 
approximately reach the nominal location of the 3:2 MMR 
(the minimum value of the period ratio is 1.508 at 170 or- 
bits). At this point in time, the planet's eccentricities increase 
to about 0.03 (inner planet) and 0.015 (outer planet). Sub- 
sequent damping of the eccentricities leads again to resonant 
repulsion, but at a much reduced pace compared to the 2:1 
MMR case shown in the uppermost row of panels in Fig- 
ure |2] The different rates of resonant repulsion are due to 
different eccentricity damping rates, which are intimately re- 
lated to the background density structure near the planets. In 
the 2: 1 MMR case, the dense gas annulus located between the 
planet gaps (see top panel in Figure |3]l provides much more 
efficient damping than in the 3:2 MMR case, where planets 
are embedded in an underdense common gap. In the latter 
case, we point out that the period ratio reaches about 1.516 
at the end of the simulation, such that the distance from the 
nominal 3:2 MMR has increased by a factor of two during the 
resonant repulsion process. This increase is consistent with 
the simultaneous factor of two decre ase in the planets eccen- 
tricities (see, e.g., lPapaloizoull201 lh . We also point out that 
no steady-state has been reached at the end of this simulation. 
For this simulation, as for all other hydrodynamical simula- 
tions of the Kepler-46 system, results are shown until the in- 
ner planet reaches r ~ 0.4, below which the proximity of the 
grid's inner edge (at r = 0.2) starts affecting the planets orbital 
evolution. 

The results of our three series of hydrodynamical simula- 
tions are summarised in Figure |4] The final period ratios are 
displayed as a function of Eo for the three different masses as- 
sumed for Kepler-46b. By final period ratio, we mean either 
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Figure 5. Time evolution of the planets' semi-major axis, period ratio, eccentricity (for the inner planet only), and resonant angles with a simple model for disk 
dispersal (see text). Disk dispersal is switched on at 1000 orbits and the evaporation timescale, T evap , is varied from 100 to 2000 orbits. Results in the upper 
panels are for Mj nncl - = 0.4Mj and Eo = 3 X 10~ 4 , and the resonant angles are for T evap = 2000 orbits. Results in the lower panels are for Mj nner = 0.6Mj and 
£o = 6 x 10~ 4 , and the resonant angles are for T cvap = 100 orbits. 



its value obtained at the end of the simulation (after between 
1500 and 2000 planet orbits), or that obtained when the inner 
planet's separation falls below r — 0.4. This figure highlights 
the three major outcomes of the simulations: (i) capture in 
2:1 MMR generally followed by rapid disk-driven resonant 
repulsion, (ii) capture in 3:2 MMR followed by slow resonant 
repulsion, and (iii) convergent migration stalling with a pe- 
riod ratio between 1 .6 and 1 .7 depending on the disk's density 
distribution within the planets' common gap. Several sim- 
ulations can actually reproduce the observed period ratio of 
Kepler-46b and Kepler-46c. The best agreement with obser- 
vations is obtained for Mj nner = 0.6Mj, our upper mass value. 
This value would imply a mean density for Kepler-46b that is 
about 20% larger than Jupiter's (that is, about 1.6 g cm -3 ). 

2.2.2. Disk dispersal model 

The results of hydrodynamical simulations in Section l2.2.1l 
show that the planets' period ratio can take a range of values 
depending on the convergent migration rate and the eccentric- 
ity damping rates during the early evolution of the system. 
The limited duration of the simulations (up to a few thousand 
orbits) and the fact that planets still migrate at the end of the 
simulations however raise the question of the long-term evo- 
lution of the system. Even when the inward migration of the 
planets gets stallecQ, further evolution of the period ratio can 
be maintained by disk-driven resonant repulsion. A steady 
state may a priori be achieved after depletion of the proto- 
planetary disk. Disk dispersal typically occurs between 10 6 
and 10 7 years after formation of the central star. Photoevap- 
oration driven by extreme ultraviolet radiation from the star 

1 This could be the case for instance when the mass of the disk inside the 
inner planet becomes smaller than the mass of the two gap-opening planets, 
a migration regime known as planet-dominated type II migration. 



opens a gap in the disk at separations of typically ~ 2 AU for 
Sun-like stars (e.g.. I Alexander & Pascuccill2012l) . For planets 
below that separation the background disk will be depleted on 
up to a local viscous timescale. Assuming a local viscous al- 
pha parameter of a few x 10~ 3 , viscous draining of the inner 
disk could operate in as short as a few thousand planet orbits. 

To assess the impact of disk evaporation on our results, we 
have restarted two simulations adopting a simple exponential 
decay of the surface density profile. This is done by solving 
<9,E = — (E — Etarget)/ Tevap m addition to the hydrodynamical 
equations, where E is the azimuthally-averaged density profile 
at restart time (from which evaporation switches on), E targ et = 
10~ 3 E is the arbitrarily small density profile in a steady state 
(taken to be non-zero for numerical convenience), and T eva p is 
the evaporation timescale. For illustrative purposes, we have 
considered three different evaporation timescales: 100, 500 
and 2000 orbits, simulations being restarted at 1000 orbits. 
Results are shown in Figure [5] The upper panels are obtained 
with Mjnnej = 0.4Mj and Eo = 3 x 10~ 4 , for which a rapid res- 
onant repulsion occurs after capture in the 2: 1 MMR. We see 
that the evolution is frozen out at a few evaporation timescales 
after restart, the planets reaching stationary eccentricities and 
semi-major axes. Note that the 2:1 resonant angles displayed 
in the fourth panel are those with T eva p = 2000 orbits. Al- 
though one of the resonant angles starts circulating from ~ 
2000 orbits, the period ratio still increases until a steady state 
is nearly reached at 3000 orbits. This shows that resonant 
repulsion requires th e libration of only one resonant angle. 
From the analysis in iPapaloizoul (1201 ll) . one can see that if 
only one of the resonant angles is effective at driving the re- 
pulsion, when it is the one containing (0\, then only e, will be 
non zero on average and slowly decreasing, being inversely 
proportional to the deviation of the ratio of the semi-major 
axes from the strictly resonant value. On the other hand, when 
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the surviving, librating resonant angle contains co , only e 
will be on average non zero. The single librating angle il- 
lustrated in the upper panel of Figure over the time inter- 
val between 2000 and 3000 orbits, for the calculation with 
fevap = 2000 orbits, involves a\ and so the eccentricity of the 
inner planet should be non zero and on average slowly de- 
creasing as is the case. Lastly, the lower panels in Figure |5] 
show the results for M; nner = 0.6Mj and Eo = 6 x 10~ 4 , for 
which convergent migration stalls at about the observed pe- 
riod ratio in the Kepler-46 system. Including disk evaporation 
stalls the migration of both planets, as expected, but it does 
not significantly affect the planets eccentricities and period ra- 
tio. The final period ratio remains very close to the observed 
value. 



2.2.3. Robustness of results 

In the course of our numerical experiments, we have tested 
the robustness of our results by varying the assumptions made 
in the physical model and by varying the grid's resolution. 
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Physical model — We performed several simulations reducing 
the disk's aspect ratio to 4% and 3% while keeping the 
same viscous alpha parameter. Although not presented here, 
these simulations show again a variety of outcomes, but 
none reproducing the observed period ratio of the Kepler-46 
system. Decreasing the disk's aspect ratio makes the gaps 
carved by the planets deeper. When convergent migration is 
rapid enough for the planets to cross their 2:1 MMR and to 
evolve in a common gap, we find that the planets generally 
lock themselves into the 3:2, 4:3 or even 7:5 MMR with very 
little departure from exact commensurability. No significant 
resonant repulsion is found to occur in these simulations as 
the density in the common gap is too low to damp the planets 
eccentricities efficiently. These simulations thus suggest 
that partial gap opening is required to obtain a period ratio 
consistent with the observations. Furthermore, as the gaps 
structure plays a prominent role in the final evolution of the 
planets' period ratio, we carried out additional simulations 
where the planets were first held on fixed circular orbits for 
500 orbits before migrating. This preliminary stage allowed 
enough time for each planet to build up a gap with a stationary 
density structure. The results of the simulations are shown 
in the upper panel of Figure [6] for Mj nner = 0.6Mj and three 
values of Eo. Starting with pre-evolved gaps changes the 
initial rate of convergent migration. Notwithstanding quite 
different initial evolutions, two out of these three models end 
up following the same evolution with or without pre-evolved 
gaps, and reach the same final period ratio at the end of the 
simulations. The model with Eo = 6 x 10~ 4 reaches a slightly 
smaller period ratio (decreased from 1 .65 to 1.6 when starting 
with pre-evolved gaps). 

Grid resolution — The effect of resolution was checked by 
doubling the number of grid cells along each direction for 
a few simulations. The results of three high-resolution runs 
(800 x 1600) are compared to the nominal-resolution runs 
(400 x 800) in the lower panel of Figure [6] Doubling the 
grid resolution slightly changes the initial rate of convergent 
migration, leading to slightly different initial evolutions. As 
for the above series of simulations with pre-evolved gaps, the 
model with capture in 2:1 MMR followed by resonant repul- 
sion is not affected by doubling the resolution. The two mod- 
els with common gap evolution seem again more SUSCepti- 
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Figure 6. Ratio of orbital periods for Mj nner = 0.6Mj and different values 
of the unperturbed surface density (Eo)- The upper panel compares our fidu- 
cial case where planets migrate from the beginning of the simulations (dotted 
curves) with the case where planets are first held on fixed circular orbits for 
500 orbits before migrating (solid curves). The lower panels shows the over- 
all good convergence of our results when increasing the number of grid cells 
from 400 x 800 (nominal resolution) to 800 x 1600 (high resolution). 



ble to a variation in the initial convergent migration, but lead 
to small differences in the final outcome. Interestingly, both 
high-resolution models with Eo = 6 and 7 x 10~ 4 yield a pe- 
riod ratio ss 1.68, in very good agreement with the observed 
value. 

2.2.4. Disk-mediated resonant dynamics 

We come back in this paragraph to the resonant repulsion 
mechanism driven by disk-planet interactions. This mecha- 
nism concerns a resonant planet pair that evolves at roughly 
constant total angular momentum, but with energy dissipa- 
tion caused by disk-driven damping of the eccentricities. Such 
configuration leads to a divergent evoluti on of the orbits an d 
thus to an increase in the period ratio ( Papal oizoul 1201 lb . 
However, it is possible that a divergent evolution of the period 
ratio caused solely by disk-planet interactions, and not by res- 
onant repulsion, could occur. To explore this possibility, we 
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Figure 7. Time evolution of the planets' semi-major axis, orbital period ratio, eccentricity and 2 : 1 resonant angles for Mj mlel = 0.6Mj and Zo = 2 x 10 . The 
panels compare the resonant dynamics in two cases: (i) the planets do not feel each other's gravity (direct interaction off) and (ii) the planets feel each other's 
gravity (direct interaction on, our fiducial case). The resonant angles are those obtained with direct interaction off, the planets resonant dynamics being mediated 
by the disk in this case. 



performed a few hydrodynamical simulations where the direct 
gravitational interaction between the planets was switched off, 
which should prevent a priori the planets from coupling res- 
onantly. The results of one of such simulations are presented 
in Figure Q and compared to a similar run where the plan- 
ets' direct interaction is included. The inner planet's mass is 
Mnner = 0.6Mj and the initial disk's density Eo = 2 x 10~ 4 
at r = 1, a configuration that leads to capture in 2:1 MMR 
followed by resonant repulsion when the direct interaction be- 
tween planets is allowed for. Surprisingly enough, we observe 
a very similar evolution when direct interaction is switched 
off: the planets period ratio increases from about 2.05 to 2.09 
between 800 and 1800 orbits. Equally surprising is the in- 
crease in the outer planet's eccentricity during the initial stage 
of convergent migration, as well as the simultaneous libra- 
tion of one of the 2:1 resonant angles (that associated with 
the periastron longitude of the outer planet and its eccentric- 
ity). These results strongly indicate that the planets are reso- 
nantly coupled although their mutual gravitational interaction 
is switched off. The reason for this is the disk: planets indi- 
rectly interact through the mediation of the disk. Each planet 
excites indeed one or two wakes (see, e.g., Figure|3j that prop- 
agate throughout the disk at the planet's angular frequency. 
The proximity of the planets implies that each planet gravita- 
tionally interacts with the other planet's wakes. When direct 
interaction is accounted for, the planets resonant dynamics in- 
volves direct planet-planet interactions and wake-planet inter- 
actions. Both types of interactions do not necessarily work the 
same way. For instance, in the simulation with no direct in- 
teraction, the inner planet's eccentricity does not change and 
the resonant angle associated with the periastron longitude of 
the inner planet circulates. 

3. APPLICATION TO KEPLER'S MULTI-PLANETARY 
SYSTEMS 

3.1. Introduction and strategy 

In Section [2] we have presented results of hydrodynamical 
simulations modeling the early evolution of planets Kepler- 
46b and Kepler-46c as they were embedded in their parent 
protoplanetary disk. The variety of period ratios obtained in 
our simulations raises the question of whether disk-planet in- 
teractions could partly explain the diversity of period ratios in 
Kepler's multi-planetary systems. 

3.1.1. Tidally-driveii resonant repulsion ? 

As highlighted in Figure Q] Kepler's multi-planetary sys- 
tems show a clear tendency for planet pairs to feature pe- 



riod ratios slightly greater than strict commensurability. For 
instance, there are almost twice as many planet pairs with 
a period ratio between 2.0 and 2.1 than between 1.9 and 
2.0. Eccentricity damping due to star-plan et tidal interactions 
could provide slow resonant repulsion (iPapaloizoul 1201 ll 
iLithwick & Wul2012tlBlitvgin & Morbidellil2013l) . Although 
the efficiency of tidal dissipation remains largely uncertain, it 
is thought that efficient tidal circularization requires the in- 
ner planet to orbit its star in less than a few days. As can be 
seen in Figure Q~]many of Kepler's candidate systems feature 
an inner planet beyond 10 days, for which tidal circularization 
is unlikely to have caused significant resonant repulsion. For 
these systems we propose that resonant repulsion driven by 
disk-planet interactions could have played a prominent role in 
determining the observed period ratios. 

3.1.2. Disk-driven type I migration ? 

The vast majority of planets in Kepler's multi-planetary 
systems have physical radii between 0.1 and 0.3 Jupiter 
radii, which places them in the super-Earth to Neptune-mass 
range. From the mass-radius diagram of Kepler candidates 
followed up by radial velocity, the median mass of a 0.2 
Jupiter-radius planet is about 10 to 20 Earth-masses (see, e.g., 
exoplanet . eu). Based on generally considered tempera- 
tures and viscosities in protoplanetary disks (corresponding to 
/; = 0.05 and a ~ a few x 10~ 3 ) such planets are expected to 
experience type I migration. Recent N-body experiments by 
iReinl (120121) have shown that convergent migration of type-I 
migrating planets yields very little departure from strict com- 
mensurability in contrast with Kepler's data. He showed that 
the inclusion of stochastic forces in addition to type I migra- 
tion could reproduce the observed distribution of period ratios 
quite nicely. Such agreement has been obtained assuming a 
fixed migration timescale and a fixed amplitude of stochastic 
forces for all Kepler's multi-planetary systems, which remains 
quite uncertain given the expected diversity of physical prop- 
erties in protoplanetary disks. Still his results indicate that 
type I migration alone may not be able to account for Kepler's 
diversity of period ratios. 

3.1.3. Disk-driven resonant repulsion of gap-opening planets? 

Based on our results of hydrodynamical simulations of the 
Kepler-46 system, we suggest that some of Kepler's planet 
candidates in multi-planetary systems could have opened par- 
tial gaps in their parent disk, thereby escaping the type I 
migration regime. This could be the case if, for instance, 
planets formed and/or migrated in regions of low turbulent 
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activity (dead zones; see, e.g., iFleming & Ston3l2003h or if 
the disk's aspect ratio takes smaller values than commonly 
adopted. To address this possibility, we carried out an hy- 
drodynamical simulation with two planets in the super-Earth 
mass range. The inner planet's mass M mnal = 15M®, the outer 
planet's mass M ou ter = 13M® and the central star is one Solar 
mass. The planet-to-primary mass ratios are therefore q\ = 
4.4 x 10~ 5 and ^ =4x 1CP 5 , that is an order-of-magnitude 
smaller than those of Kepler-46b (median value) and Kepler- 
46c. In order to get similar gap depths in the super-Earths and 
Kepler-46 simulations, we adopted values of h and a such 
that th e dimensionless p arameters in the gap-opening crite- 
rion o f lCrida et all (120061) take the same values. For the super- 
Earths simulation, this yields h = 0.023 and a ~ 2.3 x 10~ 3 . 
All other simulation parameters are otherwise identical, ex- 
cept the grid resolution that we increased to 600 x 1200. The 
unperturbed disk's surface density is Eo = 8 x 10~ 5 . The re- 
sults of this simulation are displayed in Figure [8] Planets lock 
themselves in the 3:2 MMR from about 1200 orbits. The pe- 
riod ratio reaches exactly 1.5 at about 1500 orbits and subse- 
quently goes up to 1.54 at 5000 orbits. Resonant coupling, 
fast eccentricity damping and nearly-stalled convergent mi- 
gration lead again to efficient disk-driven resonant repulsion. 
Interestingly, the lower panel in Figure [8] highlights that the 
planets do not form a common gap for periods ratios > 1.5 
in contrast to the Saturn-mass planets in the Kepler-46 sys- 
tem (compare wit h the l ower panel in Figure [3}- As already 
noticed in Section |2.2.1| the presence of the high-density gas 
annulus between the planetary gaps favors efficient resonant 
repulsion. 

3.2. Customized three-body simulations 

The results of Section [3. 1 ,3| support the idea that disk-driven 
resonant repulsion of partial gap-opening super-Earths may 
lead to period ratios that substantially differ from nominal res- 
onant values. This mechanism is intimately related to the time 
evolution of the planets eccentricities in a state of resonant 
coupling, which in turn depends on the planets masses, the 
disk's density and the time at which disk dispersal switches 
on. The survey of such a large parameter space is out of 
reach of hydrodynamical simulations, but three-body simu- 
lations with customized prescriptions for convergent migra- 
tion, eccentricity damping and disk dispersal can help illus- 
trate the diversity of period ratios that can be achieved. This 
is the strategy that we have adopted. For illustration purposes 
we have considered fixed planet masses: Mj nne r = 15M®, 
A^outer = 13M® andM* =M Q (as in Section lXOT l. 

3.2.1. Method 

Our three-body simulations solve the equations of motion 
for the two planets and the central star using a standard tech- 
nique. Planets are assumed to be coplanar. Disk-planet inter- 
actions are incorporat ed by applying ap propriate dissipative 
forces (for details, see lPapaloizoull20lll) . To model the con- 
vergent migration of partial gap-opening planets, we adopted 
a simple prescription that reproduces the results of hydrody- 
namical simulations prior to resonant repulsion. Denoting by 
R the ratio of semi-major axes (R = a uter/<2innerX we found 
good agreement using 
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Figure 8. Results of hydrodynamical simulation with two super-Earth-mass 
planets opening partial gaps in their disk. The time evolution of the planets 
semi-major axis, orbital period ratio, eccentricity and 3 : 2 resonant angles 
are shown in the upper panels. The lower panel displays the disk's surface 
density at 3000 orbital periods. 
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where pr and T are fitting parameters. The expression in 
Eq. (Q~|i tends to a constant value for ( > T, which mimics 
the fact that convergent migration is found to dramatically 
slow down or even stall before resonant repulsion sets in. In 
Eq. ©, pr denotes the period ratio at which convergent mi- 
gration approximately stalls, and T the time at which this oc- 
curs. The good agreement between this migration prescription 
and our hydrodynamical simulations is illustrated in Figure|9] 
Note that in the three-body simulations, disk-driven migration 
is applied to the outer planet only. 
Damping of the planets' eccentricity is modeled as an ex- 
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Figure 9. Time evolution of the ratio of planets semi-major axis obtained in three hydrodynamical simulations (solid curves). The planet masses and unperturbed 
surface density Eo are indicated in the top-right comer of each panel. The dashed curves display our fitting formula for the planets convergent migration prior to 
disk-driven resonant repulsion, given by Eq. iff]. 



ponential decay with a constant damping timescale to orbital 
period ratio. The simplicity of this damping prescription com- 
pared to the above migration prescription is primarily meant 
for illustration purposes, and to minimize the number of free 
parameters. We display in Figure [lO] the results of a simula- 
tion with pr = 1.50 and T = 1000 orbits, which mimics the 
convergent mi gration obtained in the hydrodynamical simula- 
tion in Section [3.1.3l (see right panel of Figure|9]l. The eccen- 
tricity damping timescale is 20 planet orbits, which is found to 
give results very similar to those of the hydrodynamical sim- 
ulation (compare with the upper panels in Figure|H}. We point 
out that the eccentricity peaks obtained near the 2:1 and 3:2 
MMRs are about a factor of 2 to 3 smaller in the three-body 
simulation. The reason for this difference is not clear. As 
shown in Section l2.2.4l it is possible that in the hydrodynami- 
cal simulations wake-planet interactions also increase the ec- 
centricities near resonances. Despite these differences, most 
importantly the evolution of the planets' period ratio com- 
pares fairly well in the three-body and the hydrodynamical 
simulations. They differ by about 25% at 5000 orbits. 

We also incorporated a simple model for disk dispersal, 
which is taken to occur between 10 4 and 10 6 orbits. Since 
planets are assumed to be already formed at the beginning of 
the simulations, this range of timescales is meant to account 
for different planet formation timescales. Once disk dispersal 
switches on, the eccentricity damping timescale is progres- 
sively increased as an exponential growth with charact eristic 
timescale T eva p = 2000 orbits. As discussed in Section 12.2.21 
this is a reasonable upper value estimate for the viscous drift 
timescale at a few tens of days (the typical present orbital pe- 
riod of Kepler's candidate planets) assuming a conventional 
viscous alpha parameter of a few times 10~ 3 . 

3.2.2. Results 

We display in FigureQT]the results of a series of three- body 
simulations using the model described in Section 13.2.11 Re- 
call thatMjnner = 15M©, M oute r = 13M ffi andM* =M©, which 
can be thought of being representative of Kepler's sample of 
multi-planetary systems (the planets' mass ratio is arbitrary, it 
takes th e sam e value as in the hydrodynamical simulation in 
Section [3.1.3b . Planets are initiated with a period ratio of 2.4. 
The prescription for convergent migration given at Eq. ([T] is 
used with pr = 1.50 and T = 1000 orbits. The value of pr 
is meant to explore the diversity of period ratios that can be 



1.8F 



1.6 



'ra 1 .4 
E 

I 

9 1 2 - 
in 



= 15 M Ea 
= 13 M F , 



1.0 



1000 2000 3000 4000 
Time (orbits) 




1000 2000 3000 4000 
Time (orbits) 











6 








M,„„,= 15 M E „ 1B 




0.008 






M„, = 13M„ - 


5 

3 


0.006 








CM 


0.004 








J? 3 

CO 

ii 2 


0.002 








3- 

1 





1000 2000 3000 4000 
Time (orbits) 



1000 2000 3000 4000 
Time (orbits) 



Figure 10. Results of a three-body simulation with prescribed convergent 
migration and eccentricity damping. The expression given by Eq. (TJ is used 
for convergent migration with X = 1000 orbits and pr = 1.5. The eccentricity 
damping timescale is 20 orbital periods. 



achieved by disk-driven resonant repulsion away from the 3:2 
MMR. The eccentricity damping timescale is varied from 20 
to 200 orbits, and the time at which disk dispersal switches 
on from 10 4 to 10 6 orbits (see Section |3. 2. \\ . Given the du- 
ration of disk-driven resonant repulsion, our results are found 
to have very little dependence on the "convergent migration 
timescale" X. Simulations were carried out over 1.2 x 10 6 
orbits, and the period ratio time-averaged over the last 1000 
orbits of the simulations is indicated by the color bar. For 
a given eccentricity damping timescale, the later disk disper- 
sal occurs, the longer planets eccentricities are damped, and 
therefore the larger the final period ratio. At a given time 
prior to disk dispersal, the shorter the eccentricity damping 
timescale, the larger the final period ratio, again as expected. 
For the shortest damping timescales that we have considered, 
period ratios can reach ~ 1.75. In this case, eccentricities are 
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Figure 11. Results of three-body simulations with prescribed convergent mi- 
gration and eccentricity damping. The central star has a Solar mass. The mass 
of the inner and outer planets is Mj nncr = 15ME ar ih and Mj mlel = 131^4^^, re- 
spectively. Convergent migration is modeled using the expression in Eq. (T} 
with pr = 1.5 and T = 10 3 orbits. The eccentricity damping timescale is var- 
ied from 20 to 200 orbits (y-axis), and it remains fixed until disk dispersal is 
switched on. The timescale at which disk dispersal switches on ranges from 
10 4 to 10 6 orbits after the be ginning of the simulations (x-axis), and the evap- 
oration timescale is set to 2000 orbits. The orbital period ratio obtained after 
1 .2 x 10 orbits, and time-averaged over the last 1000 orbits, is colour-coded. 

small enough so that planets can evolve continuously from 
the 3:2 MMR to the 2:1 MMR (that is, the 2:1 resonant an- 
gles start to librate when the 3:2 resonant angles no longer 
librate). Had we taken even shorter damping timescales or 
delayed disk dispersal even more, resonant repulsion would 
have led to period ratios exceeding 2. 

The main conclusion that can be drawn from Figure QT] is 
that partial gap-opening super-Earths may experience signif- 
icant disk-driven resonant repulsion away from the nominal 
3:2 MMR. After planets get resonantly coupled, departure 
away from commensurability is controlled by the time evo- 
lution of the planets' eccentricities, that is by the evolution of 
the disk's mass density. Under very favorable circumstances 
(particularly efficient eccentricity damping maintained over a 
particularly long timescale), period ratios can easily exceed 
nominal resonant values. Although we have considered fixed 
planet masses for illustration purposes, we believe that disk- 
driven resonant repulsion is a generic mechanism that may 
have likely occurred among Kepler's multi-planetary systems. 

4. CONCLUDING REMARKS 

The multi-planetary systems detected by the Kepler mis- 
sion have very diverse architectures, which reflect the wide 
range of conditions under which planets may form and evolve. 
The quasi-coplanar orbits in Kepler's multiple systems sug- 
gest that interactions between planets and their parent proto- 
planetary disk should play a prominent role in shaping such 
diverse architectures. In particular, close planet pairs feature 
a variety of orbital period ratios and many of them are not 
in mean-motion resonance. Convergent migration of multi- 



ple planets due to disk-planet interactions may lead to reso- 
nant systems in dis ks where turbulence is low enough not to 
disrupt resonances dKetchum et al.ll201 It IPierens et al.ll201 It 
lReinll2012h . Resonant planet pairs formed by disk-planet in- 
teractions may experience divergent evolution of their orbits 
due to tidal orbital ci rcularization (iPapaloizou & Terqueml 
120 lOt |Papaloizoull201 lb . This so-called tidal resonant repul- 
sion may increase the orbital period ratio of close-in planet 
pairs quite substantially d epending on the efficiency of star- 
plane t tidal interactions (|Papaloizoul 1201 It iLithwick & Wul 
120 12t iBatygin & Morbidellill2013l) . Tidal resonant repulsion 
could help account for the diversity of period ratios in Ke- 
pler's close-in planet pairs. 

We have shown in this paper that disk-planet interactions 
could also form planet pairs with a variety of period ratios. We 
have first focused on the Kepler-46 planetary system, in which 
the outer t wo Saturn-mass plane ts have a period ratio very 
close to 1.7 (iNesvorny et al. 2012). We have presented results 
of hydrodynamical simulations modeling the early evolution 
of planets Kepler-46b and Kepler-46c as they were embedded 
in their parent disk. For typical disk temperatures and tur- 
bulent viscosities, both planets are expected to open a partial 
gap around their orbit. A rapid convergent migration causes 
the planets to merge their gaps and to evolve into a common 
gap. Depending on the density's structure inside the com- 
mon gap, the planets' convergent migration is found to stall 
with a variety of period ratios between 1.5 and 1.7 (see Fig- 
ure 0J. The observed period ratio of Kepler-46b/c (« 1.69) 
can be reproduced without the planets being in mean-motion 
resonance. 

Our results of simulations also highlight that damping of 
the planets eccentricity by disk-planet interactions can sig- 
nificantly increase the period ratio of a resonant planet pair 
away from commensurability, a mechanism that we term disk- 
driven resonant repulsion. For instance, one of our low-mass 
disk models for Kepler-46 in which the planets get captured in 
the 2: 1 MMR may reach period ratios up to 2.3 in as short as a 
few thousand orbits. The physical process behind disk-driven 
resonant repulsion is essentially the same as in tidal resonant 
repulsion, except that the damping of the planets eccentrici- 
ties is mediated by the disk rather than by the star. While tidal 
resonant repulsion can be particularly relevant to close-in res- 
onant planet pairs, disk-driven resonant repulsion may apply 
to distant pairs. Partial gap-opening is a key ingredient for ef- 
ficient disk-driven resonant repulsion as it significantly slows 
down the convergent migration of the planets before their cap- 
ture in resonance. Smooth dispersal of the gas disk is found 
to maintain period ratios to non-resonant values. 

Many planets in Kepler's multiple systems are in the super- 
Earth to Neptune-mass range. Such planets may open partial 
gaps if they are located in a region of their parent disk where 
the local aspect ratio and/or the level of turbulence are rather 
low. In contrast to the Saturn-mass planets in the Kepler-46 
system, super-Earth-mass planets are not expected to merge 
their gaps and to evolve into a common gap, unless a very fast 
convergent migration bring them to very short mutual separa- 
tions. We have shown instead that super-Earth-mass planets 
are good candidates for efficient disk-driven resonant repul- 
sion. Depending on how fast disk-driven resonant repulsion 
operates, and when disk dispersal switches on, the orbital pe- 
riod ratio of two super-Earths may increase from 1.5 to 1.55 
or even more over typical disk lifetimes (see Figure [TT}. 

We note that a significant divergent evolution of sub-Jovian 
planets orbits was obtained in some of the hydrodynamical 
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simulations of lRein et al.ld2010l) . The similarity between their 
results of simulations and ours (convergent migration nearly 
stalled close to resonance, followed by slow divergent evolu- 
tion of the planets orbits alongside the damping of their ec- 
centricity) sugge sts that the divergent evolution obtained in 
iRein et alJ (120101) is likely to arise from disk-driven resonant 
repulsion. 

The TTV signals of Kepler's short-period planet pairs in- 
dicate that so me of these planets h ave small but finite free 
eccentricities dLithwick et al.ll2012l) . Free eccentricities get 
quickly damped by tidal orbital circularization. However, the 
presence of a disk of planetesimals or a non-smooth gas disk 
dispersal could give planets some free eccentricity. Close-in 
planet pairs that seem to have undergone resonant repulsion, 
but which have finite free eccentricities, like Kepler-23b/c and 
Kepler-28b/c, could indicate that repulsion was primarily me- 
diated by disk-planet interactions rather than by star-planet 
tidal interactions. 

We foresee at least two directions to pursue work on disk- 
driven resonant repulsion. Firstly, the effects of disk turbu- 
lence, which we have modeled with a constant viscosity, need 
to be addressed carefully. The levels of stochastic fluctuations 
required to disrupt the process as planet masses decrease, and 
to produce some residual free eccentricity should be explored. 
Secondly, a Monte-Carlo approach with three-body integra- 
tions using simple prescriptions for the effects of disk-planet 
interactions should be used to compare the distribution of or- 
bital elements among Kepler's planet pairs with the synthetic 
distributions predicted by disk-driven resonant repulsion. 



ship of the University of Cambridge. Hydrodynamical simu- 
lations were performed on the Darwin Supercomputer of the 
University of Cambridge High Performance Computing Ser- 
vice using Strategic Research Infrastructure Funding from the 
Higher Education Funding Council for England. 
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